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We study the competition between intersite and local correlations in a spinless 
two-band extended Hubbard model by taking an alternative limit of infinite dimen- 
sions. We find that the intersite density fluctuations suppress the charge Kondo 
energy scale and lead to a Fermi liquid to non-Fermi liquid transition for repul- 
sive on-site density-density interactions. In the absence of intersite interactions, this 
transition reduces to the known Kosterlitz-Thouless transition. We show that a new 
line of non-Fermi liquid fixed points replace those of the zero intersite interaction 
problem. 

PACS numbers: 71.27.+a, 71.10+x, 71.28.+d, 74.20.Mn 



Typeset Using REVTEX 



1 



The two-band extended Hubbard model is a realistic starting point both for the high T c 
cuprates and for many heavy fermion systems. The model contains a strongly correlated 
band and a weakly correlated one. At the phenomenological level, the low energy properties 
of the conventional heavy fermions (such as CeCu^ and UPt$) are well described by the 
Fermi liquid theory [Jj, while those of certain novel /—electron materials and the high T c 
cuprates appear not. The theoretical question, then, is: under what conditions do electron 
correlations lead to a non-Fermi liquid in this model? 

Recently, some progress has been made on the understanding of this model HH]. In 
the large dimensionality (D) limit, the local density- density interactions alone are found 
to cause Kosterlitz-Thouless type quantum phase transitions from a Fermi liquid to 
non-Fermi liquid metallic states. In the spinful case, the resulting non-Fermi liquids have 
spin-charge separation |||6|]. The large D limit has the advantage that local correlations 
are treated in a dynamical fashion, but the disadvantage that all intersite interactions re- 
duce to Hartree-Fock terms; no spatial fluctuations survive. For physical systems in finite 
dimensions, intersite RKKY or Superexchange type interactions are expected to compete 
with local correlations Unfortunately, for "paramagnetic" phases far away from spa- 

tial ordering transitions, there exists no controlled theoretical method to address such a 
competition. From the large D point of view, one way to recover spatial fluctuations is the 



perturbative 1/D expansion [|10j. The practicality of this procedure is unclear at this stage. 
An alternative route is a loop expansion with the requirement that the D = oo results be 
recovered at the saddle-point level, as has been constructed in models with certain forms of 



quenched disorder [[11] . For clean systems, this turns out to be difficult to formulate. 



In this paper, we take an alternative large D limit to study the interplay between local 
correlations and short range spatial fluctuations in the two band extended Hubbard model. 
We introduce an explicit intersite density-density interaction term and scale the interaction 
strength in terms of the dimensionality such that its fluctuation part survives the large D 
limit. This procedure leads to an impurity embedded in a self-consistent fermionic bath 
and a self- consistent bosonic bath. We will study primarily the spinless version of the model 



for which asymptotically exact results can be derived analytically The results are of direct 
relevance to the charge sector of the spin-charge-separated intermediate phase of the spinful 
extended Hubbard model 0,|6[|. 

The model is defined by the following Hamiltonian, 

H = YH e d n <k + Mdi + h.c.) + V:n Ci : : n di :] 

i 

+ H Mfe + KV 2 ) : n * : : n ^ W 

<ij> 

The spinless d— electrons are dispersionless, with an energy level e° d . The spinless c— elec- 
trons have a hopping matrix t,ij. < ij > labels the nearest-neighbor sites. The t— term 
describes hybridization, and the V— term the on-site density-density interaction. The in- 
terna is an intersite "charge RKKY" repulsive interaction. The standard large D limit [0] is 
taken with t <i j > scaled to be of order l/v2D and i><y> of order 1/D. In this limit, only 
the static part of the u^— term survives. Here, we scale f<y> = and take the large D 
limit keeping vq fixed. As shown below, this limit is well defined when we retain only the 
dynamical density modes. We achieve this through normal ordering, : n := n— < n >. We 
focus on states without long range ordering, for which the Hartree terms can be absorbed 
by the chemical potential. 

We first give a general formulation of this alternative large D limit. Following the 
standard procedure, we divide the Hamiltonian into local and intersite parts, 

# = + YsiKa + Ka) ( 2 ) 

i <ij> 

where hi is on-site, h t)i j = tijip\ipj, and h V)i j = (vij/2) : :: n^. :. For generality, we 
develop the formalism with the ip\ fields for possible multi- bands or spin-components; for 
the Hamiltonian ([[]), ht,ij = tijc\cj, and h V)i j = (vij/2) : :: n^. :. Within a path- 
integral representation, we divide the action for the lattice model into S = So + + AS", 
where So, S^°\ and AS" are the actions associated with ho, X^o ^« + Y,<ij^o>(ht,ij + h V; ij), 
and J2i(ht,oi + h v ,oi + H.c), respectively. Integrating out all the degrees of freedom except 
at site leads to the following effective action, S e ^ = S — E^Li( — l) n ^in^~) where 
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<>^°) denotes connected correlation functions in terms of S^°\ < AS n >)?> modifies the 
on-site action with terms that involve n site-0 density /fermion modes. The coefficients 
of these operators are retarded and are given, for each n, by n-operator correlation func- 
tions with respect to S^°\ The order in 1/D of each term can be determined through a 



cumulant expansion of the original lattice Hamiltonian We find that, to the lead- 

ing order in 1/D, no interference between h tt ij and h v>i j terms is allowed, except in lo- 
cal decorations. This absence of interference implies that, n-point correlation functions 
have the usual dependence on 1/D. As a result, for all n > 2, < AS n >( ' vanishes as 
D —>■ oo 0. The absence of interference also leads to separate Dyson equations for Xij 
and Gij in terms of their respective effective cumulants. This in turn implies that, = 
i >(''; (' j- Xij — XuXijXjji @ij = — GiiG' i0 Go G' 0j Gjj, and xif — Xij — XuXioXooX'ojXjji 
where G'y = T, P aths'tihGi 1 i 1 ti 1 i 2 G l2l2 ...Gi n i n t lnj , x'ij = Y,paths v iiiXhh v iihXhh---XUn v inj, an d 
[i, li, l 2 , In-,]} labels a non-self-retracing path from site i to site j. It then follows that, 

ij ~ V ~ ^iO^Oj/ <J"00 

xtj 1 = Xij ~ XioXoj I Xoo (3) 

From Eq. (^), <: :>(°)=< rij >' ' — < n« > is of order 1/D. Hence, in the large D limit, 
< AS >(°> also vanishes; only < AS 2 > < £ > > survives, leading to 

S eff = - [ dr / /3 c^r / ^ t (r)^ - 1 (r-r , )^(r , ) 



where 



o Jo 

-i 



+ : rty : (t)xo(t - O : : {r')\ + S (4) 



g 1 {iuj n )— — tjot^jG^ ^ 



lUJ r , 



Xo \iVn)= V i0 V 0jX\j (ifn) (5) 



IJ 

-K;„ w <o», 

y 

The Dyson equations of the lattice Hamiltonian implies that G^ is still determined by the 
on-site self-energy alone; Xij is determined by the on-site effective cumulant, or, equivalently, 
the on-site part of the vertex function irreducible in terms of both the particle hole bubble 



and the single line. All these local quantities can be calculated directly in terms of 
S e ff . Therefore, given an S e ^ , we can calculate all Gij and Xij an d hence, via. Eq. (0), 
g^iiun) and Xo fan)- The self-consistent equations (f|,|5],|3]) indeed close. They define the 
dynamical mean field equations associated with our alternative large D limit. We note 
that related mean field equations arise in the metallic spin-glass problem As usual, 

the retarded g^ 1 can be represented in terms of a self-consistent non-interacting fermionic 
bath 0. Similarly, an additional non- interacting bosonic bath can be used to represent 
the retarded Xo 1 term. We have therefore a self-consistent Anderson-like impurity model 



coupled to screening bosons [14| 



We now apply this formalism to the extended Hubbard model ([]]). The self-consistent 
impurity model is a resonant-level model with an additional bosonic bath, 

H-^l = E%dld + H + t(cJdo + h - c ) + V : n co n do : 

+ E W q p\ Pq + Y,F q :n do : (p q + p\) (6) 

where E d = e d — p, and H describes the local Cq electron coupled to a non-interacting 
fermionic bath whose dispersion, together with the parameters associated with the bosonic 
bath, W q and F q , are determined from the self-consistency eqs. (|5|). In the remainder of this 
paper, we focus on taking the large D limit with a Bethe lattice. This has the advantage 
that the bare density of states is bounded. For vq = 0, it was shown both analytically 
and numerically (T^j that the density of states of the self-consistent fermionic bath is 



non-singular in the metallic regime. We find this continues to be the case for finite Vq. This 
then allows an asymptotically exact analytic analysis of our self-consistent problem. Note 
that the spectral function of the bosonic bath is arbitrary (and is in fact non-ohmic in the 
non- Fermi liquid case, see below). 

The asymptotically exact analysis is carried out by writing the partition function of H^JJ p 
in a kink-gas representation. The procedure parallels what was detailed in Ref. |3[]. We only 
note that, the effect of the additional bosonic bath on the action of a particular history can 
be treated by a time-dependent canonical transformation, for arbitrary form of the boson 



spectral function. The resulting kink-gas action is, 

S(r 2n , n) = -2ti ln(y) + El" 1 )' h (n+i - n)/^ 

i 

+ ]T(-l)^'[2e info - n)/£ + fffa - n)] (7) 

Here, [r 2n , ri], for n = 1,2,..., labels a sequence of kink events along the imaginary 
time axis. The fugacity is given by y — tp c , where p c is the density of states of the self- 
consistent fermionic bath at the Fermi energy, e = (1/2) [1 — (l/7r)(tan~ 1 (7cp c (l — n^V) + 
tan~ 1 (irp c n ( iV))] 2 is the stiffness constant. The logarithmic interaction among the kinks is 
mediated by the fermionic bath, h = E^q is the symmetry breaking field. Finally, K(t) 
describes an additional long-range interaction of the kinks induced by the bosonic bath. Due 
to the self-consistency, K{t) is determined entirely by the full local susceptibility, 

d 2 K{r)/dr 2 = vlx{r) (8) 

Anticipating the possible l/r a form for the local susceptibility, we first solve the kink-gas 
action (^) with 

K(r)=X[(r/^r a -l]/(2-a) (9) 

for a fixed a value. We will consider the case with vanishing renormalized symmetry breaking 
field, h*. We have derived the RG equations for a close to 2, 

dyjdlni= y[l-( e + X/2)] 

de/dlni = - Aey 2 

d\/dlni= A[(2-a)-V] 

dh/dlni = h(l - 2y 2 ) (10) 

valid for small h. The RG flow diagram is given in Figure 1. The flow in the y — e plane (the 
dotted lines) describes a Kosterlitz-Thouless transition. It goes to a strong coupling fixed 
point when e < e cnt = 1, and to a line of weak coupling fixed points when e > e crii . The flow 
in the y — A plane (the dashed lines) are those of the Ising model with a long range l/r a 
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interaction [18 , 19] . There exists an unstable fixed point at (y*,X*) = (\/2 — a/2, 2), which 



describes a second order phase transition. Close to the origin, the separatrix has the form 
\ sep « 2[y/(2 — a)] 2 ~ a . The flow goes to a strong coupling fixed point for A < X sep , and to 
a weak coupling fixed point when A < X sep . In the language of the long-range Ising model, 



a = 2 is the lower critical range, and a = 3/2 the upper critical range ||18|| . The RG flows in 
between these two planes, denoted by the solid lines, interpolate between these two limits. 
There exist a line of weak coupling fixed points with y* = 0, A* = oo, and a finite e*. 

These results can be used to solve our self-consistent problem. We focus on the mixed 
valence regime only. We specify the phase diagram in terms of the three-dimensional pa- 
rameter space spanned by g t = tp , gy = [1 — (2 / 7r)ian _1 poV/2)}, and g v = poVo, where po 
is the bare density of states of the conduction electrons at the Fermi energy. The Kosterlitz- 
Thouless transition at g v — describes the charge-Kondo effect [pOfl . The critical value 
gy lt corresponds to an attractive V cnt such that p c V cnt = — {2 / Ti)tan[(^/2 — l)7r/2]. When 
gv < gy 11 -, i.e., — V < V^ T " tt , the solution is a Fermi liquid with the usual form for the local 
density susceptibility, 

X(r) ~ (1/A*r) 2 for r » 1/A* (11) 

where A* denotes the charge-Kondo energy scale, which acts as the renormalized Fermi 
energy. As V approaches V c from the Fermi liquid side, A* vanishes in the Kosterlitz- 



Thouless fashion, A* w (p c )- x exp[-l/^e crU - e]. For g v > g^ u , i.e., -V > V crit , the 
solution is a line of non-Fermi liquids with the connected local density susceptibility, 

X (r) * (1/ Pc rr (12) 

The exponent a is interaction-dependent, increasing from to 2 as one moves away from 
the critical point |T6| , |r? |. 



The intersite interaction vq modifies the phase diagram in several ways. Consider first 
9v > 9y %t ■ The line of fixed points of the v§ = problem becomes unstable. In terms 
of the parameters appropriate for the kink-gas action Eq. (^), the RG flow is towards an 



infinite value of the A coupling. Usually, one cannot specify the resulting fixed points when 
a coupling constant flows towards infinity. Remarkably, in our case we can determine the 
correlation functions in the new fixed points due to the special feature of the solution to the 
kink-gas action Eqs. (|7|, ^): the local susceptibility has an algebraic time dependence with 
an exponent identical to that of the range of K(r) pi|Ji7[] . The RG flow is towards another 
line of fixed points with an infinite A*, a finite e* , and a vanishing hybridization y*. The 
connected local susceptibility remains to have the algebraic 1/t°* form, with the exponent 
a* entirely determined by e*. 

For gv < gy lt , we are able to establish the existence of a phase transition as vq is 
increased. First, the solution must be a non- Fermi liquid for sufficiently strong v Q . Suppose 
that the solution is a Fermi liquid, with a renormalized Fermi energy E*. Then the local 
susceptibility has a long time (1/E*t) 2 dependence for r longer than 1/E*. At times up 
to 1/E*, if the local susceptibility decays slower than 1/t 2 , then the corresponding K{r) in 
this intermediate time range has the form Eq. (|9]) with a < 2 and A = v 2 A, where A is 
the prefactor in the decay of the local susceptibility. The scaling equations ( |T0D imply that, 
as long as e + A/2 > 1, the fugacity decreases in this time range. If the local susceptibility 
decays faster than 1/t 2 , one can still choose a v sufficiently large such that the fugacity does 
not increase up to the time scale 1/E*. In both cases, at times beyond 1/E*, the kink-gas 
action corresponds to a Coulomb gas with a stiffness constant e' = e(£ = 1/ E*) + (v / '\/2E*) 2 
and a small value of fugacity y'. When vq becomes sufficiently large, (e', y') will lie in the 
weak coupling side of the RG flow, making E* — 0. Therefore, the Fermi liquid solution 
cannot occur. The only self-consistent solution in this large vq regime is one similar to what 
we found for the gv > gy %t case, characterized by a y* = 0, a finite e* , and an infinite A*. 
The local susceptibility has the form of Eq. (|I2]), with the exponent being smaller than 
two. Physically, the intersite density-density interactions provide charge-screening, which 
contribute to the orthogonality effect [§,0 . In the mixed valence regime, this orthogonality 
helps realize the weak coupling fixed point with incoherent charge excitations. 

For sufficiently small vq, on the other hand, the Fermi liquid solution is stable. This 
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can be seen by a vq— expansion around the vq = solution. We replace K(t) in the kink- 
gas action Eq. (|7|) by what we would get from Eq. (|j) with the local susceptibility x( r ) 
of the v = problem, Eq. ([EL]). At the time scale £' = 1/A*, e(£') < e(£ ) < 1, and 
K{t) = (vq/A*) 2 ln(r/^'). From Eqs. @, the kink-gas action at £ > £' is a Coulomb gas 
with a stiffness constant e' = e(£') + (t> /v^2A*) 2 . For sufficiently small vo, e' < e cnt = 1. 
In this range, the Fermi liquid solution is self-consistent. Self-consistency, however, will 
modify the renormalized Fermi energy, making it unlikely that the phase transition is of the 
Kosterlitz-Thouless type. The precise nature of the transition is beyond the reach of our 
RG formalism. The schematic phase diagram is given in Fig. 2. 

As a result, non- Fermi liquids with self-similar correlation functions occur even for re- 
pulsive values of the on-site density-density interaction. 

In the non-Fermi liquid case, that the renormalized A* is infinite is one indication that 
the ground state cannot be the "paramagnetic" metallic state. The fact that the Ising 
model corresponding to the kink-gas action Eqs. has a divergent free energy at zero 

temperature for a < 1 implies the same physics. The precise nature of the ordering depends 
on the details of the band structure and the intersite interactions. Our results applies at 
temperatures above the ordering temperature. 

The extension of the self-consistent equations and the kink gas action to the spinful 
extended Hubbard model is straightforward. The form of the scaling equations implies that 
our results carry over to the charge sector of the spin-charge-separated intermediate phase. 

After the completion of this work, we learned of the independent work of Kajueter and 
Kotliar [22| who constructed related mean field equations in the context of a spinless one- 
band fermion model with semi-circular density of states and found no numerical evidence for 
non-Fermi liquids in that model. We thank V. Dobrosavljevic, A. Georges, and A. Sengupta 
for stimulating discussions. Q.S. was supported by an A. P. Sloan Fellowship, and by NSF 
Grant No. PHY94-07194 at ITP, UCSB. 
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FIGURES 

FIG. 1. The RG flow of the kink-gas action Eqs. (g, |) for a <~ 2 and vanishing renormalized 
symmetry breaking field. 

FIG. 2. The phase diagram of the Hamiltonian Eq. ([!]) in the mixed valence regime. The 
dashed line is schematic. 
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